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Abstract: We have demonstrated recently that RHIC heavy ion data concerning particle produc- 
tion at small and intermediate transverse momenta become very transparent, if one realizes that 
there are two sources of particles production: the corona - due to the interactions of the periph- 
eral nucleons of either nucleus, and the core - representing the high density central zone. We 
extend our analysis to SPS heavy ion data. Again, we find that the nontrivial centrality depen- 
dence of particle production (being stronger than at RHIC) is simply due to an increasing corona 
contribution with decreasing centrality. The core contribution shows no centrality dependence 
(concerning particle ratios), and can be described in exactly the same way as at RHIC, with one 
exception: there is somewhat less collective radial flow at SPS compared to RHIC. Apart of this 
(and a trivial volume effect) , the core portions produced in different colliding systems at different 
energies (RHIC, SPS) look the same. 



One of the spectacular results from SPS 
heavy ion research is the very strong enhance- 
ment of strange (in particular multi-strange) 
baryons compared to proton-proton scattering. 
Although this goes qualitatively into the right 
direction, having in mind the production of a 
quark-gluon-plasma |1|, a quantitative under- 
standing is still lacking. 

We do not really claim to explain these data, 
we simply want to stress that they show an 
almost trivial behavior, if one takes into ac- 
count the "corona effect": the peripheral nu- 
cleons of either nucleus essentially perform in- 
dependent pp or pA-like interactions, with a 
very different particle production compared to 
the high density central part. For certain ob- 
servables, this "background" contribution com- 
pletely spoils the "signal", and to properly in- 
terpret the data, we need to subtract this back- 
ground. 

In order to get quantitative results, we em- 
ploy exactly the same method as described re- 
cently 1 2 1 for analyzing RHIC heavy ion data: we 
use EPOS 1 3 1, which has proven to work very 
well for pp and dAu collisions at RHIC and for 
pp scattering at SPS. EPOS is a parton model, 



so in case of a nuclear collision there are many 
binary interactions, creating partons, which 
then hadronize employing a phenomenological 
procedure called string fragmentation. Here, 
we modify the procedure: we have a look at 
the situation at an early proper time t , long 
before the hadrons are formed: we distinguish 
between string segments in dense areas (more 
than p a segments per unit proper volume), from 
those in low density areas. We refer to high 
density areas as core, and to low density ar- 
eas as corona. In order to make a quantitative 
statement, we adopt the following strategy: the 
low density part will be treated using the usual 
EPOS particle production which has proven to 
be very successful in pp and dAu scattering 
(the peripheral interactions are essentially pp 
or pA scatterings). For the high density part, 
we simply try to parameterize particle produc- 
tion, in the most simple way possible. To do so, 
we consider the core contributions sepatately 
in different longitudinal segments. Connected 
core regions in a given segment are considered 
to be clusters, whose energy and flavor content 
are complete determined by the correspond- 
ing string segments. Clusters are then consid- 
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ered to be collectively expanding. We assume 
that the clusters hadronize at some given en- 
ergy density Ehadr. having acquired at that mo- 
ment a collective radial flow, with a linear radial 
rapidity profile from inside to outside, char- 
acterized by the maximal radial rapidity y ra( j. 
In addition, we impose an azimuthal asymme- 
try proportional to the initial space eccentricity. 
Hadronization then occurs according to covari- 
ant phase space. For more details and param- 
eters see [2|. 

So we employ exactly the same procedure 
as we did for RHIC, with even the same pa- 
rameters - up to two exceptions, concerning 
the parameters r and y Ta ^. Whereas at RHIC 
the final results are insensitive to variations of 
r (in the range l-2fm), this is no longer the 
case at SPS, the reason being the finite reac- 
tion time T rea (time it takes for the two nuclei to 
pass through each other), which is somewhat 
more than a fermi. So we use r = r rca , rep- 
resenting the minimum possible value. Since 
we have smaller initial core densities at SPS 
compared to RHIC, we also expect a smaller 
radial flow, so we take the freedom to use 
//,,„[ as a free parameter, fixed by comparing 
to SPS data (the only free parameter when 
going from RHIC to SPS!). We actually use 
y rad = 0.60 at SPS, instead of 0.83 at RHIC. 

In the following, we will discuss results for 
PbPb collisions at 158 GeV. In fig. ^ we com- 
pare the core contribution corresponding to 
a central (0-5%) PbPb collision (which means 
purely statistical hadronization, with flow) with 
pp scattering. We plot m t spectra of pions, 
kaons, protons, and lambdas, the nuclear 
spectra are divided by the number of binary 
collisions (according to Glauber). Apart of the 
shape differences, the most striking feature is 
the fact that the yields for the different pp con- 
tributions are much wider spread than the core 
contributions, even more than at RHIC. In pp 
(string fragmentation), lambdas (and even more 
the multi-strange baryons) are very much sup- 
pressed compared to pion production, and this 
suppression is even stronger at SPS compared 




mt-m 

Figure 1: Invariant yields l/2nm t dn/dydm t of 
pions (red), kaons (blue), protons (green), and 
lambdas (yellow), for the core contribution 
corresponding of a central (0-5%) PbPb colli- 
sion (full lines) and proton-proton scattering 
(dashed). The PbPb spectra are divided by the 
number of collisions. 

to RHIC (since the string masses are smaller). 
In fig. |21 we compare the full contribution (core 
and corona) of a central (0-5%) PbPb collision 
and proton-proton scattering . The PbPb pion 
line crosses the pp line at large p t , which is due 
to the fact that the core contains particles from 
pA like collisions, where we have a "Cronin- 
enhancement" due to parton ladder splitting, 
discussed in detail in (3), which has been in- 
troduced in order to understand dAu scatter- 
ing at RHIC. In fig. |3j we plot the relative 
contribution of the core (relative to the com- 
plete spectrum, core + corona) as a function 
of m ( - m, for different particle species. For 
central collisions, the core contribution domi- 
nates largely (around 90%), whereas for semi- 
central collisions (40-50%) the core contribu- 
tion decreases, giving more and more space for 
the corona part. The precise m t dependence 
of the relative weight of core versus corona de- 
pends on the particle type. 

We are now going to study SPS PbPb data 
(158 GeV). In fig. HI we plot the centrality 
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Figure 2: Invariant yields l/2irm t dn/ 'dydm t of 
pions (red), kaons (blue), protons (green), and 
lambdas (yellow), for the full contribution (core 
+ corona) of a central (0-5%) PbPb collision (full 
lines) and proton-proton scattering (dashed). 
The PbPb spectra are divided by the number 
of collisions. 



dependence of the particle yield per partici- 
pant (per unit of rapidity), for different parti- 
cle species, the data |4] |H] together with the 
full calculation (upper diagram), and the full 
calculation compared to the core contribution 
(lower diagram). The complete calculation fol- 
lows quite closely the data. Whereas central 
collisions are clearly core dominated, the core 
contributes less and less with decreasing cen- 
trality. In fig. |SJ we consider the correspond- 
ing particle ratios, as a function of centrality, 
for the core contributions. We show the ratios 
of different particles, with respect to 7r _ or K s . 
The rations are practically flat, apart of some 
decrease for very small participant numbers. 

So our first important conclusion: af- 
ter subtracting the "corona background", 
the interesting part, the core contribution, 
shows an extremely simple behavior: there 
is no centrality dependence, the systems 
are simply changing in size. The partici- 
pant number is certainly not a good measure 
of the volume of the core part, this is why 
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Figure 3: The relative contribution of the core 
(core/(core+corona)) as a function of the trans- 
verse mass for different centralities (0-5%: red, 
40-50%: blue. Upper figure: pions (full) and 
kaons (dashed). Lower figure: protons (full) and 
lambdas (dashed). 



the overall multiplicities per participant de- 
crease with decreasing centrality. This is 
the same conclusion as in \2\ for AuAu at 
RHIC. And not only the SPS core contri- 
bution is as simple as the RHIC one, it is 
even parameterized with the same parame- 
ters, apart of somewhat more flow at RHIC! 

What makes the measured centrality depen- 
dence look complicated is the mix of core and 
corona, depending on the particle species. Ac- 
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Figure 4: Upper diagram : Particle yields per 
participant in PbPb scattering at 158 GeV as 
a function of the number of participants, for 
(from top to bottom): ir~ , K + , K~, K s , A, A, 
S, 5, ft, Cl. We show data (points) (H| to- 
gether with the full calculation (core + corona, 
full line). The three upper -most curves refer 
to 47rresults, all others to central rapidities. 
The leftmost points are pp calculations (col- 
ored points) and pp or pBe data (black sym- 
bols). Lower diagram : full calculations, same 
as in upper plot (full lines) compared to the core 
contributions (dotted). 



tually the Q are the simplest, since here the 
corona contribution is negligible. Here, all 
the participants which contribute to the corona 
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Figure 5: Particle ratios as a function of cen- 
trality, from top to bottom: K + /tt~, K~/tt~, 
k/K s xO.\ , A/K s xO.l, E/K s xO.l, 3/K„ xO.l, 
Q./K s xO.l, Cl/K s xO.l. The full lines are the 
calculations for the core contribution alone, the 
dashed lines are horizontal lines, just to guide 
the eye. 

(being more and more frequent with decreasing 
centrality) do not contribute at all to the fl pro- 
duction. This is why the yields drop so strongly 
at small participant numbers. This also reflects 
the fact that the core volume is not proportional 
to the number of participants. 

To demonstrate consistency, we check in 
the following m t spectra for different hadron 
species, for PbPb collisions at 158 GeV, see figs. 
|6] and [TJ The calculations represent the full 
contribution, core plus corona, whose relative 
contributions (at least for kaons and lambdas) 
can be obtained from fig. |3j So it is actually a 
non-trivial superposition of these two contribu- 
tions, hiding the shape of the interesting part, 
the core. 

It should be noted that NA49 |8| reports 
much less lambda production at central rapid- 
ity than NA57, and correspondingly our calcu- 
lation of dn/dy(y = 0) (being close to the NA57 
values) is about 40% above the NA49 results 
(whereas the transverse mass spectra almost 
agree!). 
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Figure 6: Transverse mass spectra of K s (upper 
plot), A (middle), A (lower) in PbPb scattering at 
different centralities. From top to bottom: 0- 
5%, 5-11%, 11-23%, 23-40%, 40-53%. Lines 
are full calculations, points are data |6|. The 
different curves are displaced by factors of 10. 

In fig. |8j we show rapidity and transverse 
mass spectra of different hadron species for 
central CC and SiSi at 158 GeV. Again we just 
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Figure 7: Transverse mass spectra of S (upper 
plot), 5 (middle), fi+fj (lower) in PbPb scattering 
at different centralities. From top to bottom: 0- 
5%, 5-11%, 11-23%, 23-40%, 40-53%. Lines 
are full calculations, points are data The 
different curves are displaced by factors of 10. 

show the full calculations (core plus corona) 
together with the data, and we find a similar 
good agreement than in case of PbPb (using 
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the same parameters!), with the exception of 
the cf> transverse momentum spectrum, which 
is somewhat harder in the calculations. Con- 
cerning the relative weights of core and corona, 
CC and SiSi is very similar to peripheral PbPb, 
in the sense that the corona is relatively more 
important than in central PbPb (the surface ef- 
fect is bigger in small nuclei than in big ones) . 



To summarize: we have discussed the in- 
fluence of the corona contribution (occurring 
in the periphery of nuclear collisions) in PbPb, 
SiSi, and CC collisions at the SPS. We provide 
a realistic treatment of the corona, by using 
a model which works excellently for pp and 
pA. We can provide a parametrization of the 
core, such that the complete calculation (core 
+ corona) provides a good fit of the observed 
particle spectra in heavy ion collisions at the 
SPS. Our core results represent to some extent 
"background corrected results", removing the 
undesired corona contributions from the spec- 
tra. The core shows an extremely simple 
centrality and species dependence: actually 
there is no centrality dependence, all par- 
ticle ratios are constant. In other words, 
the cores in all these different systems (CC, 
SiSi, PbPb at different centralities) are all 
identical, apart of a trivial volume differ- 
ence. Even more remarkable: the SPS re- 
sults (concerning core) are identical to the 
RHIC results, with one exception: there is 
30% more radial flow at RHIC. So what re- 
ally makes the experimental results compli- 
cated is the non-trivial admixture of corona 
contribution to the core. 



All the results concerning SPS in this paper 
have been obtained with exactly the same pro- 
cedure (and same parameters) as in a previous 
study of RHIC results 2|, with just one excep- 
tion: the radial flow parameter. All the central- 
ity and system size dependence is purely deter- 
mined by geometry. 
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Figure 8: Rapidity spectra of identified hadrons 
in central SiSi (upper, left) and central CC (up- 
per right) and transverse mass spectra in cen- 
tral CC (middle) and central SiSi (bottom). In 
each figure, from top to bottom: tt + , tt~ , K + , 
K~ , A, <j). Lines are full calculations, points are 
data. 
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